Scalar-isoscalar states in the large-iV c Regge approach 
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Scalar-isoscalar states (J PC = ++ ) are investigated within the large- V c Regge approach. We 
elaborate on the consequences of including the lightest /o(600) scalar-isoscalar state into such an 
analysis, where the position of /o(600) fits very well into the pattern of the radial Regge trajectory. 
Furthermore, we point out that the pion and nucleon spin-0 gravitational form factors, recently 
measured on the lattice, provide valuable information on the low-mass spectrum of the scalar- 
isoscalar states on the basis of the scalar-meson dominance in the spin-0 channel. Through the fits 
to these data we find m CT = 450 — 600 MeV. We compare the predictions of various fits and methods. 
An analysis of the QCD condensates in the two-point correlators provides further constraints on the 
parameters of the scalar-isoscalar sector. We find that a simple two-state model suggests a meson 
nature of /o(600), and a glueball nature of /o(980), which naturally explains the ratios of various 
coupling constants. Finally, we note that the fine-tuned condition of the vanishing dimension-2 
condensate in the Regge approach with infinitely many scalar-isoscalar states yields a reasonable 
value for the mass of the lighest glueball state. 
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I. INTRODUCTION 

The history and status of the cr-meson has been quite 
vacillating (for reviews see e.g. [ll-Q and references 
therein). A scalar-isoscalar state with a mass of ~ 
500 MeV was originally proposed in the fifties .4] as an 
ingredient of the nucleon-nucleon force providing satura- 
tion and binding in nuclei. Along the years, there has 
always been some arbitrariness in the "effective" or "fic- 
titious" cr meson mass and the coupling constant to the 
nucleon, partly due to the lack of other sources of infor- 
mation. For instance, in the very successful Charge De- 
pendent (CD) Bonn NN-potential any partial wave 
2S+1 L,/-channel is fitted with a different scalar-isoscalar 
meson mass and coupling. The cr-meson was also intro- 
duced as the chiral partner of the pion to account for 
spontaneous breaking of the chiral symmetry Q. The 
lack of confidence in its existence motivated taking its 
mass to infinity, yielding the non-linear sigma model [7|, 
which is the modern starting point for the Chiral Pertur- 
bation Theory Q. 

During the last decade, the situation has steadily 
changed, and the cr-meson has been finally resur- 
rected [9| , culminating with the inclusion of the ++ res- 
onance in the Particle Data Group review (PDG) [T3| as 
the /o(600) state, seen as a 7T7t resonance. It has wide- 
spread values for the mass, 400 — 1200 MeV, and for the 



width, 600 - 1200 MeV [ll|. A rigorous definition of the 
cr as a 7T7T resonance requires that it be a pole of the 
7T7T scattering amplitude in the (J, T) = (0,0) channel in 
the second Riemann sheet in the Mandelstam variable s. 
Within such a framework, the uncertainties have recently 
been narrowly sharpened with a benchmark determina- 
tion based on the Roy equations with constraints from 
the chiral symmetry [12| , yielding the value 1 



tT ct /2 = 441 



+ 16 



z272Zi2 MeV. 



(1) 



The analysis of Ref. [13] yields a result with somewhat 
higher 77v, 473 ± 6 ± 11 - «257 ± 5 ± 2 MeV (with the 
errors statistical and systematic, respectively), while the 
unitarized Chiral Perturbation Theory (yPT) give s a bit 
lower value of the mass, 4011^-7277^26 MeV [Til. Nev- 
ertheless, various determinations of the pole mass agree 
with the values ([1]) within the uncertainties. These ac- 
curate determinations make somewhat tricky the orig- 
inal question on what a mass should be used a priori 
within a meson-exchange picture, due to the very large 
width of the resonance. Moreover, the determinations 
mentioned above do not imply necessarily the standard 
assignment of the linear sigma-model where one takes 
(cr, 7?) as chiral partners in the (5,5) representation of 
the chiral SU (2) R £g> SU (2) L group. A priori, the cr state 
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1 We use this definition for the pole mass in the yfs variable. A 
better one is s CT = — iTo-nia-, which coincides with the previ- 
ous one in the narrow resonance limit. 
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could also belong to the (0,0) representation. Admit- 
tedly, the heated debate on the nature of the a meson 
is not completely over, with the tetraquark [l5[ (see also 
[HI) or glueball interpretations [ItJ considered (see, e.g, 
Ref. [l[ and numerous references therein). 

Our goal in the present paper is to point out that the 
large- N c Regge models may provide valuable insight into 
this problem. Clearly, the numbers in Eq. (JIJ represent 
the values for N c = 3. We recall that in the large- N c limit 
QCD can be mapped onto a theory of infinitely many 
stable and non- interacting mesons and glueballs [la , Il9j . 
In the standard large-A c counting, meson and glueball 
masses are both 0(N°) whereas their widths are C(iV ( r 1 ) 
and (D(N~ 2 ), respectively. However, there is also a large- 
r-suppressed mass shift which makes a literal use of 
Eq. (P) somewhat questionable. Actually, lattice calcu- 
lations have shown that QCDqo is generally not too far 
from QCD at N c = 3 (see, e.g., Ref. [13] and references 
therein), at least for certain observables. The lightest 
scalar-isoscalar state (commonly denoted by a), due to its 
unusually large width, seems to be an exceptional case, 
and in this paper we try to address its nature and phc- 
nomenological consequences within the framework of the 
large- N c Regge models. 

Large- Ac-motivated investigations for the irir or NN 
interactions keep only the lowest ++ state, and hence 
constitute the low-energy analyses (see e.g. Ref. [2l[ for 
an early investigation) . The studies of the 7T7t scattering, 
based on a large- A c scaling of the A c = 3 parameters and 
unitarized via the Inverse Amplitude Method (IAM) ap- 
plied to the ;\TT amplitudes, provide an A c dependence 
of resonances which, regarding the a state, depends on 
the details of the scheme used. While the one-loop cou- 
pled channel approach [22j yields a very wide range of vn a 
and a large width (in an apparent contradiction with the 
standard large- N c counting [18|, [l9( when extrapolated to 
N c >• 10), the two-loop approach [23| produces a large 
shift (by a factor of 2) of m a when going from N c = 3 to 
N c = oo (the corresponding mass shift is small in the case 
of the p meson) . One should note the large uncertainties 
of the two-loop IAM, documented in Ref. 24j. Al ong 



similar large-A c counting, the Af-matrix approach [25] 
gives m a ~ l^prifi^ = 0(\fW c ). The Bethe-Salpeter ap- 
proach [26| yields an estimate m a ~ 500 MeV [27| at 
N c — > oo. A more refined analysis using the large- N c 
consistency conditions between the unitarization and res- 

„ = 0{N~ 1 ) ' 



14] 



onance saturation suggests m p — rr 

The large- N r . analysis of the NN force was carried out 
in Refs. [28|, [29| at the quark level and in Ref. [30] at 
the hadronic level, providing a justification of the meson- 
exchange picture. Moreover, it indicates what mesons 
are the leading ones, namely the infinite tower of a, tt, 
p, lo, Ai, etc., states. Meson widths enter as relative 
1/N 2 corrections to the potential, on equal footing with 
many other effects (the spin-orbit force, relativistic dy- 
namics, or the inclusion of other mesons), independently 
on how large the a width is in the real N c — 3 world. 
In this regard the large-A c analysis may provide a re- 



liable determination of the asymptotic value of the a 
mass at N c — > oo. A fit to the 1 5'o leading- N c poten- 
tial yields m a = 501(20) MeV [13]. Moreover, it was 
shown in Ref. [31] within the large-A c framework, that 
the correlated 2n exchange corresponds to a Yukawa po- 
tential where the corresponding Yukawa mass differs by 
a 1 /N c correction to the pole mass or equivalently to the 
Breit-Wigner mass (the difference between both masses is 
0(N~ 2 ) [32( ). Finally, we should also note that using an 
alternative large- N c counting [33|, [34| , the scalar meson 
behaves as a qqqq state |35j and its width is (D(N~ 2 ). 

In the present work we take into account an infinite 
number of scalar-isoscalar states. In order to gather in- 
formation on their coupling to hadrons and the vacuum, 
we use a Regge model for the radial trajectories and study 
correlation functions between the energy-momentum ten- 
sor. Surprisingly, we find in Section [IT] that quite natu- 
rally all the scalar-isoscalar states can be described by 
a single radial Regge trajectory with half the standard 
slope. The mass of the a state can be deduced from this 
trajectory as the mass of the lowest state. As a conse- 
quence, there seems to be no obvious difference between 
mesons and glueballs, as far as the spectrum is concerned. 
Despite this surprising ordering with a sort of a radial 
quantum number, we should warn the reader that states 
of different nature may not necessarily be characterized 
by a single quantum number and other interpretations 
are considered in the literature. For example, /o(980) 
is often regarded as qqqq, /o(1370) as mainly uu + dd, 
/o(1500) as mainly a glueball, and /o(1710) dominantly 
as ss [l|. In such a case, there is no trajectory reveal- 
ing the sigma meson mass, and the observed regularity 
would be accidental. 

Existing lattice studies of hadronic matrix elements of 
the energy-momentum tensor, the so-called gravitational 
form factors of the pion and nucleon, provide a very 
valuable complementary information on the a properties, 
with the large- N c interpretation behind (Section II VI) . 
Unfortunately, the data are too noisy as to pin down 
the coupling of the excited scalar-isoscalar states to the 
energy-momentum tensor. Nevertheless, useful informa- 
tion confirming the mass estimates for the cr-meson can 
be extracted (Sect. IIV[) . It is observed from the lattice 
data that the a mass grows with the value of the light 
current quark mass in a natural way. 

It is interesting to discern the nature of the a state 
from an analysis of a truncated spectrum. The minimum 
number of states, allowed by certain sum rules and low 
energy theorems, is just two. In Section [V] we undertake 
such an analysis, which suggests that /o(600) is a qq me- 
son, while /o(980) is a glueball. Finally, in Section |VT1 we 
address the rather fundamental issue of the existence of 
dimension-2 condensates within the context of the Op- 
erator Product Expansion of QCD. Such objects appear 
naturally in the Regge approach and could not vanish if 
only a finite number of states were kept. Actually, we 
show that the dimension-2 condensate may vanish when 
infinitely many states are considered, in which case a 
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FIG. 1: (Color online) Radial Regge trajectory corresponding 
to the squared mass of all J PC 
listed in the PDG tables [H| 



scalar-isoscalar states 
The four heaviest ++ states 



are not yet well established and are omitted from the PDG 
summary table. The error bars correspond to the errors in 
the determination of the square of mass, listed in [3(|. The 
straight line is the result of Fit I with half the standard Regge 
slope, according to Eq. with a and m a from Eq. Q. The 
band reflects the uncertainty in the parameters of Eq. @. 



fine-tuning condition on the a mass is needed. Finally, 
in Section [VIII we come to our conclusions. 



II. SCALAR-ISOSCALAR REGGE 
TRAJECTORIES 



In this section the mass of the cr-meson is determined 
by extrapolating the upper part of the J PC = ++ spec- 
trum in the large- N c Regge approach. In Ref. [37j the 
analysis of the radial and rotational Regge trajectories 
was carried out. In Ref. [38[ the scalar sector was stud- 
ied in a greater detail, including the modeling of the me- 
son widths. Two parallel radial trajectories could then 
be identified, including three states per trajectory, which 
has resulted in rather poor statistics. In Fig. [1] we show 
the squared mass of all nine J PC = ++ isoscalar states 
listed in the PDG tables (36[, given for completeness in 
Table HI Note that the four heaviest ++ states are omit- 
ted from the PDG summary table and in principle need 
confirmation. As we can see, an approximate straight line 
behavior looks quite appealing. To improve the statistics, 
we propose to fit with half the standard slope, thus in- 
cluding all the PDG scalar-isoscalar ++ states together 
in the mass formula 



M s {nf 



a 

2 U 



(2) 



where a = 2na, and a is the string tension. The same 
strategy was recently applied in Ref. [39| . 

In the first method, called Fit I, we minimize the x 2 



E 



M 



M s {n) 



AM 



(3) 



We get, after excluding the lowest n — state and using 
the 8 remaining PDG states, the central parameter values 
a = 1.35 GeV 2 and m a = 530 MeV, with x 2 /DOF = 
X 2 /(8 - 2) = 9.9. 

A non-trivial issue, not addressed in the analyses of 
Refs. [U Hl|, is the determination of errors of the pa- 
rameters in the Regge fits. Actually, the large value of 
X 2 /DOF ~ 10 3> 1 found above prevents a reliable er- 
ror determination of the model parameters a and rn a 
and calls for corrections from different sources. A com- 
mon way to overcome this difficulty is to rescale the 
weights in the \ 2 variable. In our case we simply re- 
place AM/,„ -> 3AM/, n , such that x 2 /DOF ~ 1. In this 
case we get 



1.35(5) GeV 2 



Din 



530(95) MeV (Fit I). 



(4) 



In the usual approaches (see, e.g., Refs. [38|, |40j) a dif- 
ference between the bare and physical pole masses, based 
on the iGmatrix or the Breit-Wigner formula, is made. 
This generally introduces some model dependence and 
thus there may be systematic uncertainties in the masses 
related to the precise definition of the resonance param- 
eters. 

As we can see from Eq. ([4]), the fit is quite accept- 
able and implies the value of the string tension y/a = 
463(9) MeV. We note that the phenomenological value 
is y/a = 420 MeV from a global fit to the excited meson 
spectrum (37j . 

Formula ^ is equivalent to two parallel radial Regge 
trajectories with the standard slope 



Ms-(n) 2 = an + m 2 , 
M s ,+ {n) 2 = an + ml + -. 

From here we get the mass-splitting formula 



and, in particular, 



2(m 



/o(980) 



27TCT, 



(5) 



(6) 



(7) 



which works extremely well for the value m / = 980 MeV 
and m<7 = 530 MeV, yielding 1.36 GeV 2 for the l.h.s. vs. 
1.35 GeV 2 for the r.h.s. 

Given the fact that all states fit quite naturally into 
the pattern, we see no obvious way how glueball states 
could be singled out solely on the basis of belonging to 
one of the trajectories ([5]). Also note that a glueball 
corresponds to a bound state in gluodynamics, or equiv- 
alently QCD for infinitely heavy quarks (zero active fla- 
vors), while the states obtained here correspond to the 
limit of light quarks (two or three active flavors) . 
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Resonance 


M [MeV] 


T [MeV] 


n 


M (Fit I) 


M (Fit II) 


/o(600) 


400 - 1200 


500 - 1000 





530 


556 


/o(980) 


980(10) 


70(30) 


1 


976 


983 


/o(1370) 


1350(150) 


400(100) 


2 


1275 


1274 


/o(1500) 


1505(6) 


109(7) 


3 


1516 


1510 


/o(1710) 


1724(7) 


137(8) 


4 


1724 


1714 


/o(2020) 


1992(16) 


442(60) 


5 


1909 


1896 


/o(2100) 


2103(8) 


209(19) 


6 


2078 


2062 


/o(2200) 


2189(13) 


238(50) 


7 


2234 


2215 


/o(2330) 


2321(30) 


223(30) 


8 


2380 


2359 



TABLE I: PDG values of resonance parameters [36j, com- 
pared to the fits to the radial Regge spectrum of the ++ 
scalar- isoscalar states, M% = \an + m 2 . In fits I and II the 
uncertainties are taken as the error in the mass, or as the one 
half of the resonance width, respectively. In both cases the 
lowest n = state, corresponding to the /o(600) resonance, 
is excluded from the fit. 
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FIG. 2: The A\ 2 = 2.3 and 4.6 contours, corresponding to 
the 68% and 90% confidence levels in the m CT — a plane, for 
Fit I (solid lines) and Fit II (dashed lined). The optimum 
values are indicated with the dots, filled for Fit I, and open 
for Fit II. 



As is well known [l8|, [l!|, at large iV c mesons and 
glueballs are stable states with fixed mass, Ms, n ~ , 
and small decay widths Ts, n ~ 1/N C and ~ re ~ 
spectively. Moreover, a general feature of the observed 
mesonic spectrum is that while the masses grow with n, 
the widths stay constant. Therefore, one would naively 
expect the large- N c approximation to work better for the 
higher masses of the excited states. One might incorpo- 
rate the widths by using a particular model for the spec- 
tral density, such as the Breit-Wigner parametrization. 
Rather than using such a parametrization for the finite 
width, it may be better to take the widths themselves as 
an uncertainty on the value of the mass up to 0(1/N C ). 
This is clearly an upper bound on the uncertainty of the 
resonance position. Actually, the PDG resonance val- 



ues are often quoted as the Breit-Wigner values, which, 
depending on the assumed background, changes from 
process to process and in general leads to some model 
dependence. However, a rigorous description of reso- 
nances as quantum mechanical states requires to deter- 
mine them as poles of the meson-meson scattering ampli- 
tudes in the second Riemann sheet. Then their position 
in the complex plane becomes model-independent. For 
the lowest ++ state the shift between the Breit-Wigner 
and the pole values is rather large, m a = 800 MeV vs. 
rn a = 444 MeV, so the difference is compatible with the 
corresponding half-width. In view of this discussion we 
fit Eq. ([2]) by minimizing 



E (— 



M s {n) 



W 2 



(8) 



where we have chosen the inverse of the half-width 
squared as the weight. The result of this procedure, 
called Fit II, is 



a = 1.31(12) GeV 2 , 
X 2 /DOF = 0.12 



= 556(127) MeV, (9) 
(Fit II). 



The correlation plot for the a — m a parameters is pre- 
sented in Fig. [2] for fits I (solid lines) and II (dashed 
lines). We notice the compatibility of the two methods. 

We remark that the large- iV c -motivated analysis of 
the NN-scattering in the ^-channel [13, S3 yields a 
quite similar value for the lowest scalar mass, m a = 
500(20) MeV. 

Although our fitting strategy was rather crude, the re- 
sulting Regge-like spectrum of the scalar-isoscalar states 
turns out to be very reasonable, given the simplicity of 
the mass formula of Eq. @. The crucial aspect of the 
foregoing analysis is that the /o(600) resonance corre- 
sponds to the lightest state in the radially-excited spec- 
trum of the Regge-like scalar-isoscalar states. This is 
an important point, since scalar glueballs [HI, |43[ and 
scalar-isoscalar mesons [3 [45[ have been studied within 
the AdS/CFT framework, however neglecting the possi- 
ble role of the lightest scalar in any of these approaches. 
The soft- wall version has been reviewed in Ref. [Ill |4tJ . 

In Sect. II Vl we will provide a mounting evidence on the 
coupling of the lightest state to the gravitational form 
factor of both the pion and the nucleon. This is relevant, 
as it shows that the a state should be explicitly consid- 
ered in the evaluation of the correlation functions, where 
it plays an essential dynamical role. 



III. FORMALISM FOR THE EFFECTIVE 
SCALAR FIELDS 

In this Section we review the basic formalism used in 
the following parts of the paper. A particular attention is 
paid to the A^-scaling of coupling constants of the meson 
and glueball scalar-isoscalar states. 
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A. Interpolating fields 



B. Two-point correlations 



The first non-trivial question is to look for the inter- 
polating QCD operator which triggers the J PC = ++ 
isoscalar states from the vacuum. On the lattice this has 
been done in a variety of ways , m ainly involving quark 
operators (such as qq, see also (48|). A clear candidate, 
often used for the gluonium, is given by the trace of the 
energy-momentum tensor, Q^" , which satisfies the trace 
anomaly equation (49[, 



d>*D u 



e 



(10) 



~2c7 & G ^ 



9 



m q [1 + 7m («)] qq- 



Here /3(a) = [x 2 da/d[i 2 denotes the beta function, a is 
the running coupling constant, J m (a) — d log m/d log /i 2 
is the anomalous dimension of the current quark mass 
m, the symbol denotes the dilatation current, and 
Cr" is the field strength tensor of the gluon field. The 
operator O, besides having the J PC = ++ quantum 
numbers, is renorm-invariant. In addition, in the chiral 
limit of massless quarks, m q —> 0, it is SUr(2) ® SUl(2) 
chirally invariant. This is different from the usual qq 
interpolating field, which is not chirally invariant and 
mixes under the chiral transformations with the qij^fq 
operator. Actually, under the operation q 75^, the 
operator O is chirally even while qq is chirally odd. 
In perturbation theory one has 



7m (a) 



—a 

a 
4tt 



f- 

V47T 



(11) 



where 



A) 
ft 



3 

34 Ar2 
3 c 



-N 



-N f N r 
3 ' 



N f 
N, 



— , (12) 



and A/y denotes the number of active flavors. We recall 
that in the large- N c limit a ~ 1/N C and G 2 ~ (TV 2 - 1), 
i.e., is proportional to the number of gluons. Hence, in 
Eq. 1[1 Op there is a contribution scaling as N 2 and flavor- 
independent, as well as another contribution, subleading 
in N c and scaling as N c Nf (the higher orders in, for in- 
stance, the MS scheme, generate the N 2 terms which are 
dependent on the renormalization scheme). In the chiral 
limit we only have the gluonic operator contribution to 
Eq. (|10|) . however, still some information on the quark 
degrees of freedom remains via the /3-function. Thus, in 
the chiral limit we may distinguish between the gluonic 
and quark contributions by the large- N c and Nf scaling 
behavior, 



e 



N? 



e 



(13) 



For a scalar-isoscalar particle \n) we have a non- 
vanishing matrix element to the vacuum for the dilatation 
current, 



and hence 



\D»\n)=iq»f n , 



(0|cFZ?» = (0|e|n)=m 2 /„ 



(14) 



(15) 



for the on-shell particles. From now on the mass of 
the scalar-isoscalar states is denoted as m„. The two- 
point correlation function of the energy-momentum ten- 
sor reads 



ee 



(?) 



i|T{6^(a;)e Q/3 (0)} |0). 



(16) 



The energy momentum tensor is conserved, therefore 



(Olefin) = i/„(ff"V - <f<f), 



(17) 



where the factor of 1/3 complies with Eq. f|15[>. After 
inserting a complete set of states we find 



H 



ee 



with 



q a q^U ee (q), 
(18) 



U ee (q) = i d i xe^ x (0\T{Q(x)Q(0)}\0) 



+ c.t., 



(19) 



where c.t. stands for the countertcrms. According to the 
N c counting rules, the leading part of nee scales as N 2 , 
and the next-to- leading part as N e Nf. 

The sum in Eq. ([T9"l) runs over all J PC = ++ isoscalar 
states, i.e. both mesons (m) and glueballs (g). The dif- 
ference is, however, in the scaling of the coupling con- 
stants /„ with N c . If all states are contributing and 
no cancellations occur, then the glueball must have a 
larger coupling, f 2 ~ N 2 , whereas the qq- meson has 

fn~N c N f . 



fn 
fn ' 



(glueball), 
(meson). 



(20) 



The Op erator Product Expansion (OPE) in QCD 
yields 



H(g 2 )=<? 4 



C log q 2 



(21) 



6 



where 



C. Effective Lagrangeans 



Co = - 



(22) 



In the conventional treatment insisting on the presence of 
the local gauge-invariant operators only, the sum over n 
starts with n = 2, while the admission of the dimension-2 
condensate includes the n = 1 term as well. In Sect. IVII 
we will analyze this issue further. 

We now match the QCD result (f2"Tj) to the large- N c for- 
mula (|19p . Firstly, we notice that to generate the log q 2 
term we need infinitely many states to contribute. Sec- 
ondly, approximating the scries with an integral, we find 
via the Euler-Maclaurin formula that 



f n /(dm n /dn) — > const. 



(23) 



For the Regge spectrum in Eq. @ this specifically means 
that f 2 does not asymptotically depend on n. 

The asymptotic q log q 2 behavior allows one to write 
down a twice-subtracted dispersion relation, 



U( q 2 ) =n(0)+n'(0)q 2 



IT 



dt ImlT(t) 

t 2 t — q 2 — ie 



In our case, the spectral density is 



1 



ImII(s) 



f n m n S{m n - s), 



hence we can write 



mt — q 



(24) 



(25) 



(26) 



We have imposed the condition to have a well behaved 
object at high ener gies . With this prescription the low- 
energy theorem [5~i1 l55l |56| , is written as 

Ilee(O) = £ f n ml = -4(9) = -16B V , (27) 

n 

where B v is the energy density of the vacuum. The first 
derivative with respect to q 2 yields 



n^ (O) = V/ n 2 



d 4 ra 2 (6(x)e(0)), (28) 



a dimension-2 object which cannot be expressed as an 
expectation value of a local operator. The second deriva- 
tive is equal to 



ee 



(o)=e4 

^— ' mi 



(29) 



a dimensionless object which corresponds to the gravi- 
ton wave-function renormalization and which may be di- 
vergent. In general, if we use the fact that /„ becomes 
asymptotically rt-independent and m 2 ~ n, as suggested 
by the radial Regge trajectories, the sums (|2"7|) and ((55J) 
need regularization (see Sect. IVII for an application of 
the ^-function) and the naively positive combination of 
Eq. (|28l) may in fact turn into a vanishing or even nega- 
tive contribution. 



The results of the previous subsections can be directly 
translated into the language of the effective Lagrangeans. 
In particular, we may write the PCDC (Partial Conser- 
vation of the Dilaton Current) equation, 



^(*) = £/„3>„(*) + e% 2 ), 



(30) 



and hence for the on-shell particles we may use the equa- 
tions of motion, (—d^d^ — m^)<^„ = 0, to obtain 



^/ n m 2 ^„(x)+0(^ 2 ). 



(31) 



For the off-shell particles we need to introduce the im- 
proved energy momentum tensor .57], which has a better 
high energy behavior, but the final result is the same. 

The above properties may be described with a suitable 
extension of the effective dilatonic Lagrangeans [58l - [6l| 
of the form 



n 

Here, the potential 



2i$ 4 
P2 n 



log 



fn 



(32) 



(33) 



fulfills the conditions that the minimum is at ($„) = f n 
and the mass is V"(f n ) = m 2 n . The physical field ip n is 
defined by $„ = /„ + cp n . The value at the minimum is 



(34) 



which corresponds to the vacuum energy density (|2~T1) . 

The coupling to the pion can be considered by working 
in the non-linear representation and with the interaction 
Lagrangean 



4 4" /1 



rm-K ,2 



$ 2 (^[/t^C7), 



(35) 



where U (x) — e lT '^^ 7 ' is the non-linearly transforming 
pion field and (.) represents the trace in the isospin space. 
In the vacuum, $„ = f n , we get 



4 0) - ^£<w/»<^W>> 



which leads to the condition 



(36) 



(37) 
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Since ~ N c and /„■ ~ y^ci it follows (under the 
assumption that all states arc contributing in the large- 
N c limit and no cancellations occur) that gmrnfn ~ W c °, 
or, due to Eq. pi)) . 
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■ 1/JV C (glueball), 
(meson). 



(38) 



A construction of the chiral Lagrangean in the linear re- 
alization is presented in Appendix |Bj 

The scalar fields $ n are invariant under chiral trans- 
formations. 2 Lagrangean (l3"5j) is nothing else but agen- 
eralization of the dilaton Lagrangean of Refs. |Hj, l60j | for 
the case of many dilatons. We treat here all scalars on 
equal footing, as suggested by the fact that they all fit 
into the same radial Regge trajectory. 

The tree-level decay width of the scalar n into two 
pions is given by the expression 



F($„ 7T7T) = F n7r7r 



(39) 



32tt 



1 



4m 2 



3.9 



2 ni 3 



32tt 



Because of the scaling 

r„ w ~ 1/A 2 
- l/iV c 



TLTTTZ 

r 

- 1 - ■t'l.TTTT 



(glueball) , 
(meson), 



(40) 



which displays the weaker coupling of the glueballs to the 
pion compared to the mesons. 

The Lagrangean describing the interaction of the nu- 
cleon with the scalar-isoscalar particles is 



C 



N 



N 



(41) 



with u 2 = U 5 = e^^/f* . The nucleon field N trans- 
forms non-linearly, but we can return to the linear re- 
alization by undoing the chiral transformation, = 
uN, which effectively replaces the nucleon mass Mjv by 
the combination J2 n g n NN$n- This yields the general- 
ized Goldberger-Treiman relation for the scalar-isoscalar 
channel, 



M 



N 



(42) 



derived by Carruthers in the early seventies |62J. Again, 
under the assumption that all states contribute to (|4"2")l 
and no cancellations occur, we find 



9nNN 
9nNN 



A c ° (glueball), 
V^Vc (meson) . 



(43) 



Lagrangean (14T1) fulfills also the standard Goldberger- 
Treiman relation, g^Af/v = f^g^NN- In addition, it is 
scale-invariant, such that as a consequence that is 
the Noether current generating dilatations we have the 
identity [HMJ 
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E 

n 

4f2 ^n' 



[47(*„) - <S>nV\<$> n ) 



(44) 



Upon subtraction of the vacuum part and using that 
= fn + (fin it follows that 



(N\e(0)\N) = - Y,m 2 Jn{N\ Vn \N) 



fngn 



NN, 



(45) 



in agreement with (|4"2"1) . In other words, the pole mass 
and the energy-momentum mass are the same, as re- 
quired by consistency. 

The contribution of the scalar-isoscalar degrees of free- 
dom to the ttN scattering is given by the formula 



T: b N = 6 ab u(p')u(p) 



•sr~^ Q 9mrTr9nNN 

„ ml-q 2 



(46) 



which vanishes at the threshold q = 0, in agreement with 
the Weinberg- Tomozawa result. According to Eqs. (|38|) 
and (|4"0")) . the meson contribution scales as N® and the 
glueball contribution as 1/N C , hence is relatively sup- 
pressed. 

Finally, the nucleon-nucleon central potential stem- 
ming from the exchange of the J PC = ++ isoscalar 
particles reads in the large- iV c limit 



Vb(|q|) 



n 2 

„ 1 2 + ™»' 



(47) 



Passing to the coordinate space we obtain the standard 
Yukawa potential, 

. „2 -m„r 

Vc(r) = -T 9 -^ e - (48) 



The large-iV c counting rules (f4*U|) yield Vc ~ N c from 
the meson exchange, while the glueball exchange pro- 
duces a subleading contribution ~ N®. Other tenso- 
rial components of the NN potential scale as ~ N c or 
- l/N c llil. 

The N c scaling of various quantities is collected for 
convenience in Table [TTJ 



IV. GRAVITATIONAL FORM FACTORS 



2 This is unlike the standard assignment of the linear sigma-model 
where one takes (cr, 7?) as chiral partners in the (1/2, 1/2) rep- 
resentation of the chiral SU(2)fi Cg) SU(2)l group. See also Ap- 
pendix [B] 



In this section we use the lattice data for the gravita- 
tional form factors of the pion and nucleon to place fur- 
ther bounds on the mass of the a meson. Our analysis 
is based on the fact that the spin-0 channel is dominated 
by the scalar-isoscalar mesons. 
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TABLE II: iV c -scalirig of various quantities. 



quantity 


giueball 


qq meson 


m„ 


1 


1 


/» 


N c 


vm 


r 


l/Nl 


l/No 


Qnnn 


1/N C 




g?iNN 


1 





A. Pion 

The gravitational form factors of the pion are defined 
via the matrix element of the energy-momentum tensor 
between the pion states, 

(ttV) I ©^(O) I 7r°(p)> = 

l _ 5 ab _ q^e^ + ipvp'Q^ > (49) 

where P = ^(p' + p), q — p 1 — p, and a, b are the isospin 
indices. The trace part of this form factor corresponds 
to the coupling of the spin-0, while the traceless part to 
spin- 2 states. In the present analysis we are interested in 
the former. The spin-0 gravitational form factor is given 
by 



^(p')|e(o)| 7 r a (p)) = ( 5 ab e.(^) 



(50) 



where 



e.(<Z 2 ) = lq 2 Ql(q 2 ) + - 9 2 )e 2 (g 2 ). (51) 

Due to chiral symmetry constraints, this form factor sat- 
isfies pSlL l55l. [Hal . 



A 0.52 

H 

V 

0.51 



0.50 - 




0.42 0.44 0.46 

m a (GeV) 



FIG. 3: The A\ 2 = 2.3 and 4.6 contours, corresponding to 
the 68% and 90% confidence levels, in the m a — {x) u +d plane. 
The fit is obtained from the lattice data for the pion spin-0 
gravitational form factor, as described in the text. 



exactly as Eq. (|37|) . Since according to Table |U we have 
the scaling fn^n-ww ~ , both mesons and glueballs may 
contribute to relation (j5"7| on equal footing. 
At large momenta the form factor behaves as 



',! k 7i Win f n 



If we require the constant term to vanish, we get 



^ fnQrnr-n'Wln ^* 



..(58) 



(59) 



Q 7r (q 2 ) = q 2 + 0(p 4 ). 
Two low energy theorems follow: 



e w (o) 
euo) 



0, 
1. 



(52) 



(53) 
(54) 



These conditions can be achieved by assuming a deriva- 
tive coupling 



(7r°(p')|* n |7r 6 ( P )) =S' 



ab 9n 



9 9 ' 

m l n - q z 



written as 



and thus 



e.(<z 2 ) = £ 



£ 2 2 

m£ - q z 



(55) 



(56) 



(57) 



Clearly, cancellations are necessary and we need at least 
two contributing states to satisfy the sum rule ([5^]) . This 
high-energy relation plays a similar role in the scale- 
invariance considerations to the second Weinberg sum 
rule in the chiral framework. 

The only determination of the pion gravitational form 
factor to date is the lattice calculation of Brommel et 
al. [H, H3. The data go up to t = -4 GeV 2 , however 
with substantial errors. The quark contribution to the 
energy momentum tensor, considered in (63l. Ifiij. is 



0f (x) = \q{x) 



7 M iD v + Y iD 1 - 



q(x), (60) 



from which a decomposition similar to Eq. (|49[) and the 
corresponding spin-0 component Eq. (|51[) follow. The 
data points are shown in Fig. [4] Following the standard 
notation for the moments of the pion GPD's [65| we in- 
troduce 



1 



A 20 (t) = -el(t), A 22 (t) 



em, 



(61) 
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where the symbols denote the quark parts of the grav- 
itational form factors of Eq. (J3U). A monopole fit has 
been undertaken, yielding ra a — 0.89(27) (9) MeV and 
A 22 (Q) = -0.076(5) (cf. Table C.8 on p. 105 of Ref. (H). 
By using several parameterizations it has been noted that 
the precise form of the function cannot be pinned down 
from the data unambiguously due to large uncertainties 
in A22{t). In the following we will use the low-energy 
theorem 



^22(0) =-~A m (P), 



(62) 



which provides a relatively accurate fixing of A22 at the 
origin and greatly helps the regression analysis. 

In terms of the form factors related to the Generalized 
Parton Distributions (GPD) of Ref. [H, HI], the quark 
contribution to the gravitational form factor reads 

0«(f) - -AtAUt) (63) 

Due to the multiplicative QCD evolution one has 



(64) 



where the (valence) quark momentum fraction depends 
on the renormalization scale /i. Its leading-order pertur- 
bative evolution reads 



R 



> 9 (M) 



7l 0) /(2/3o) 



(65) 



where the anomalous dimension is 7^ - > /(2/3 ) = 32/81 
for Np = N c = 3. The QCD running coupling constant 
is equal to 



a(/z) 



4?r 



A,log(/x 2 /A 2 



(66) 



QCD J 



where we take A QC d = 226 MeV for N c = N f = 3. 

Phcnomcnologically, it is known from the Durham 
group analysis [661 ]. based mainly on the E615 Drell- 
Yan data [67| and the model assumption that the sea 
quarks carry 10 — 20% fraction of the momentum, that 
(x)l +d = 0.47(2) at the scaleu = 2 GeV. The analy- 
sis of the Dortmund group [68| . based on the assump- 
tion that the momentum fraction carried by the valence 
quarks in the pion coincides with that in the nucleon, 
yields (x)^ +d = 0.4 at n = 2 GeV. The lattice data at 
the lattice spacing a\ at = O.lfm as well as a recent chiral 
quark model calculation [69J support this view. 

We recall that the large- N c analysis implies sums of 
monopoles in the form factors. The largest available mo- 
mentum transfer, t = —4 GeV 2 , obtained in Ref. f&ij ]. 
suggests that some information on the contribution of the 
excited states might be extracted. Therefore, following 
the approach already used for the electromagnetic pion 
form factor in Ref. [7(| (see also [7l| ) , we have attempted 
a Regge-like fit, 



e„(t) = tf b {t), 



(67) 



including infinitely many states, of the form 



h(t) 



(68) 



with B(x,y) — T(x)T(y) /T(x + y) denoting the Euler 
Beta function. The function (|6"5)l fulfills the normaliza- 
tion condition 



MO) = 1. 



(69) 



For x 3> y one has B(x,y) ~ ^(y)x y , hence in the 
asymptotic region of M 2 - f > (ii - l)a we find 



hit) 



a/2 
M 2 — t 



b-l 



(70) 



The result for a = 1.31 is 



0.52(3), m 



495lf|g MeV, b 



2 OA+ 1 ' 56 
Z ' Z ^-0.55- 



(71) 



As we can see, the result is fully compatible with the 
monopole (b = 2) and at present the large errors in the 
lattice data wash out any insight from the excited scalar 
spectrum, despite the large momenta 3 . 

Thus, we restrict ourselves to a simple monopole fit 



— ( x )u+d 2 



mi - 1 



(72) 



which yields x 2 /DOF ~ 2.4. Actually the large \ 2 i s 
due to incompatible values of nearby points. In such 
a situation, in order to obtain reliable estimate of the 
model parameters, we rescale the errors by a factor of 
1.5 to make them mutually compatible. Moreover, we 
enforce the low energy theorem, Eq. (|62p as a constraint 
- a possibility not directly considered in Ref. [H, l&jj . 
As a result we get % 2 /DOF ~ 1 (after the mentioned 
rescaling of the data errors) and the optimum values 

( x )l+d = 0.52(2) m a = 445(32) MeV. (73) 

In Fig. [3]we present the corresponding correlation ellipse. 
The gravitational form factor Oo(t) at the optimum val- 
ues of the parameters ([75)) is presented in Fig. S) 



B. Nucleon 

The scalar component of the nucleon gravitational 
form factor is 



(N(p')\e(o)\N( P )) = u( P ')u(p)e N (q 2 



(74) 



3 We note that a similar fit [7l| to the vector form factor using 
a = 1.2(1) GeV 2 yields m p = 775(15) MeV and b = 2.14(7), 
which can be distinguished from a simple monopole fit at the 
two-standard-deviation level. 
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-t [GeV~] 



-A 20 (0)/4 



FIG. 4: (Color online) Spin-0 gravitational form factor of the 
pion, A22{t), from the lattice calculation of Refs. [(H, [64[ ex- 
trapolated to the physical pion mass (the points with the large 
error bars), together with the monopole fit with parameters 
(|73|l . indicated by the line. The band indicates the uncer- 
tainty in the model parameters. The low energy theorem 
fixes the value of the A22 form factor at t = 0. The lower 
data points show the data for — A20/4. 



where u(p) and u(p') are nucleon Dirac spinors. In the 
meson-dominance approximation it can be written as 



NN- 



fnm 2 n 



(75) 



The nucleon mass is given by Eq. (|42l) . 

The quark contributions to the GPD's of the nucleon 
have been determined on the lattice [72|, |73| for masses as 
small as twice the physical pion mass. The decomposition 
corresponding to the energy-momentum tensor reads 



(p'|e*» = u( P >) 



2 



4Mjv 
- g^A 2 



M N 



u( P ), (76) 



where a pv = §[7^,7"] (the Bjorken-Drcll notation), the 
scalar functions are moments of the GPD's, the momen- 
tum transfer is denoted as A = p' — p, and the average 
nucleon momentum is P = (p' + p)/2. Taking the trace 
and applying the Gordon identity, 2MNu(p')^ tJj u(p) — 
u(p')(i<T^ p A p + 2P M )w(p), as well as the Dirac equation, 
(j> - M N )u{p) = and u{p')(p" - M N ) = 0, we obtain 
the following expression for the spin-0 gravitational form 
factor of the nucleon: 



M 



N 



AUt) 



t 



B q w {t) 



3f 



M 



2 C|o(*) 

N 



.(77) 



Due to the multiplicative character of the QCD evolution 
for the conserved energy-momentum tensor operator, one 
has 

£ e%(t) = "£{x)»e N (t) = (x) u+d Q N (t). (78) 




-t [GeV 2 ] 

FIG. 5: (Color online) Spin-0 gravitational form factor of 
the nucleon, Ge(t), obtained from the lattice simulations of 
Refs. [731 at the pion masses m n = 352 and 356 MeV, together 
with the monopole fit with mass m CT = 624 MeV (solid line). 
The band indicates the uncertainty in m„ of 78 MeV obtained 
with the % 2 method. 




0.4 0.5 0.6 0.7 
m<T [GeV] 



FIG. 6: The A\ 2 = 2.3 and 4.6 contours, corresponding to 
the 68% and 90% confidence levels in the m a — c plane. The fit 
is obtained from the nucleon spin-0 gravitational form factor 
obtained from the LHPC data [73J. See the text for details. 



The combination (|77p can be directly constructed from 
the results of the LHPC Collaboration where the 
data points for the Aiq, -B201 an d C20 form factors are 
provided for six values the pion mass: = 757, 681, 
595, 495, 356, and 352 MeV. The values of the corre- 
sponding nucleon masses are Mn = 1565, 1489, 1379, 
1292, 1216, and 1158 MeV. We add the statistical errors 
in quadrature. The result for the two lowest pion masses 
is shown in Fig.[SJ We note that while the A 2 o(t) function 
is determined to a better than 10% accuracy, the large 
errors in Q q N originate from the less accurate B20 (t) and 
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C2o(i) form factors. 

Similarly to the pion case we have tried to deduce some 
information from the high-energy region, t ~ —1.2 GeV 2 , 
by attempting the Regge form 

Bib — 1 m2 ~* ) 
Q N (t)=M N '~gT . (79) 

Unfortunately, again the large errors in the data do not 
allow to fix the b parameter unambiguously, and it is 
compatible with 6 = 2. Thus we stick to the simple 
monopole form, 

Qj f {t)=M N {x)^ d -^-. (80) 
m£ — t 

Since most of the lattice data are accumulated at pion 
masses significantly above the physical value, extrapola- 
tion to the physical point is necessary. At low values of 
the pion mass and t the chiral corrections are expected to 
be significant. 4 Not entering these intricacies, we assume, 
following many lattice studies, the simple dependence of 
m a on m n : 

mKm-n) = ml + c(ml - m 2 phys ) . (81) 

Hence m a (without an argument) denotes the value at 
the physical point. Similarly, 

( x )u+d( m *) = ( x )u+d + d(ml- ml phys ) . (82) 

We use simultaneously the form factors (jTT)) for the six 
pion masses from LHPC [73]. There are several ways one 
can perform the fit, all leading to similar but not identical 
results. Our procedure is as follows: We first normalize 
the form factors to unity at t = 0, by dividing a data set 
at a particular by the value of the point at t = 0. 
That way we get rid of the dependence on (x)^ +d and 
can carry out a two-parameter fit, with the parameters 
of Eq. flSIJl. The x 2 method yields 

m a = 550±2oo MeV ' c = °- 95 -a75> ( 83 ) 

with x 2 /DOF = 5.2/(87 - 2). The resulting correlation 
plot is shown in Fig. [6] 

One may perform separately the fit to the normalized 
form factors at each value of . An example for the two 
lowest pion masses, 356 and 352 McV, is shown in Fig. [5] 
The band indicates the uncertainty in m a obtained with 
the x 2 method. The result of the procedure at all values 
of is shown in Fig. [71 together with the best fit curve 
with parameters ([55)1 . The point at the physical value of 



4 Experience with the vector form factor shows that the low mo- 
mentum region is quite sensitive to chiral corrections, whereas 
the intermediate energy region is better described by the large- 
N c dynamics. 




S 0.4 - 



0.2 - 

0.0 L 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 ' 

0.0 0.2 0.4 0.6 0.8 

m K [GeV] 

FIG. 7: (Color online) Dependence of the a mass on in- 
ferred from the LHPC data [7l|. The line shows the best fit 
of Eq. H83j> . The point at physical indicates the result of 
the extrapolation. 




S 0.4 - 
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0.0 0.2 0.4 0.6 0.8 1.0 

m n [GeV] 

FIG. 8: (Color online) Dependence of the a mass on in- 
ferred from the QCDSF data [7^. The line indicates the best 
fit of Eq. (85]). 



mjr indicates the resulting extrapolation for the physical 
value of m a , Eq. (155)) . 

Taking the values of the unnormalizcd form factors at 
t = we may compute the momentum fraction carried 
by the quarks as a function of m T . The result of the \ 2 
fit is 

(1)^ = 0.447(14), d= -0.13(3). (84) 

In the case of the lattice calculations of the QCDSF 
Collaboration [72| , five points have been provided at val- 
ues ranging between < — t < 3 GeV 2 for the pion 
masses m w = 930, 760, and 550 MeV, then extrapolated 
to the physical pion mass. The corresponding nucleon 
masses are M N = 1724, 1519, and 1308 MeV. The A, B, 
and C form factors have been parametrized in Ref. [72[ 
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by a dipole form with the mass of M — 1.1(2) GeV. 5 

Note that in the presence of a heavy mass of a tensor 
meson the data for the spin-0 gravitational form factor 
become noisy at large values of the momentum trans- 
fer t due to the relative smallness of the signal and the 
propagation of errors in the formula for 9^, Eq. (|77|). 
For this reason we proceed as follows: with the provided 
dipole fits we reconstruct the data points for the spin-0 
gravitational form factor using Eq. (ITT)) at the quoted 
values of t, adding the statistical errors in quadrature. 
The resulting Q 9 N (t) becomes negative within errors at 
—t > 1.5 GeV 2 , precluding accurate fitting with the 
monopole form. Thus we determine the mass from the 
slope of the linearly-extrapolated form factor at t = 0, ac- 
cording to the formula 1/to 2 — ^Q q N (0)/Q q N (0). Figure 
[8] shows the outcome of this procedure for the available 
pion masses as well as for the data set resulting from an 
extrapolation to the physical point made in Ref. [72[ • 
The result of the fit of the form (f8Tj) yields 

m CT = 600t?g MeV, c = 0.8(2), (85) 

or mt = 0.35 GeV 2 + 0.8m 2 . These values are consistent 
within the uncertainties with the formula of the LHPC 
collaboration, to 2 = 0.387 GeV 2 + 0.91m 2 .. 

Similarly, fitting to the form factors at the origin, we 
obtain 

(x)% +d = 0.55(3), d= 0.07(4). (86) 
C. Summary of the a mass 



n 1 1 1 i 1 1 1 i 1 1 1 i 1 1 1 i 1 1 1 i 1 1 1 r 
Particle Data Group 

spectrum, Regge Fit I 

spectrum, Regge Fit II 

n gff, Regge (QCDSF data) 

k gff, monopole (QCDSF data) 
i — • — i 

N gff, monopole (LHCP data) 

slope of N gff, monopole (QCDSF data) 

200 400 600 800 1000 1200 

mo- [MeV] 

FIG. 9: (Color online) The collection of various results for 
the mass of the sigma state. The vertical band indicates the 
"benchmark" calculation of Ref. [12|]. The fits leading to the 
six bottom points are described in the text. 



The various estimates of the mass of the lightest scalar- 
isoscalar state, ttv, are summarized in Fig. [S] We note 
that the values are compatible within the uncertainties. 
We remark that the definition of the resonance mass via 
the pole position on the second Riemann sheet of the 
scattering amplitude or via the Breit-Wigner function 
may differ significantly, even by a factor of two. How- 
ever, the difference occurs at the level 1/iV 2 (32^ . 

The previous results are based on extrapolation to 
physical pion masses, a subject which has received much 
attention in recent years, as chiral logs are expected to 
play a role. Comparison of Figs. [7] and [8] provides a co- 
herent picture that the scalar mass grows with the cur- 
rent quark mass. Actually, by analyzing the heavy pion 
limit, a distinction between glueballs and mesons could 
be made. Indeed, glueballs are existing states in pure glu- 
odynamics which corresponds to the limit of extremely 
heavy quarks. The lightest ++ glueball in that limit has 
a mass of M ++ = 1.7GeV, and this would correspond to 



5 While in Ref. [72l | it is claimed that data disfavor a monopole 
shape, it is also said that the results support the assumption of 
the tensor meson dominance of the monopole form for the spin-2 
form factor. 



a physical state whose mass becomes a fixed number for 
infinitely heavy quarks. Clearly, states which become qq 
states should decouple since they have masses scaling as 
2m q at large m q . The data are just too scarce and noisy 
to make a conclusive statement in that regard. From this 
viewpoint, lattice data in the heavy-pion regime would be 
most welcome. 



V. CALCULATIONS WITH A FINITE NUMBER 
OF STATES 

In this Section we return to the derived sum rules in an 
attempt to asses the nature of the scalar-isoscalar states. 
The sum rules (|34I42I57I59I) are: 

^/nflW = 1, (87) 

n 

^ ^ fn9mr7T rn n ^, 
n 

fn9nNN = Mn, 

71 

- _L f 2 2 _ 

_ 16 / j Jn m 7i ~ € v 
n 
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FIG. 10: (Color online) Ratios of coupling constants of the 
a and /o(980) in the two-state model. All ratios are of the 
order of 1/ \/N c , indicating the interpretation of the the a as 
a qq meson and /o(980) as a glueball, cf. Table ITT1 

In addition, the dimension-2 object of Eq. (|21I28I) is 

T,fn = C 2 . (88) 

It is clear that the satisfaction of the second sum rule (|8"7) 
with a finite number of states requires the presence of at 
least two scalar-isoscalar states. The presently-available 
lattice data do not allow to extend the analysis of Sect. lIVI 
to account for two or more states. Nevertheless, one can 
perform a qualitative analysis based solely on Eq. (|87|) . 

When we saturate the sum rules (|87[) with two states, 
we get 

1 = 9a-jnvfa + 9finrff, (89) 
= ga^faml + gfTnrffmj, (90) 
M N = g a NNfcr+9fNNff: (91) 

ev = ~{flml + f f m)), (92) 
C 2 = fl + f}- (93) 

Let us identify a with /o(600), the state / with /o(980), 
and treat m a as a variable. From the decay widths 
= 540(20) MeV I2j] and T fo ^ = 70(30) MeV 
we infer (up to the sign) with the help of Eq. (|3"9"|) the 
m^-dependent values of the coupling constants g a -mv and 
gfirir- Then, equations (|89| are used to find the constants 
fa and //. We choose the convention where both of these 
constants are positive. Next, assuming g a NN = 9.5(5), 
we obtain gjNN from Eq. (|9"Tjl . Finally, from Eq. (|9"2")l 
we compute the vacuum energy density. The results are 
displayed in Figs. [TOJ HQ and[T2J 

The most important qualitative feature can be in- 
ferred from the analysis of the ratios of coupling con- 
stants. Figure H0l shows the ratios of the constants fa/ ff, 
9fmr/g<7-Kni and g/NN I g&NN ■ All these ratios are of the 
order \/^/W c ~ 0.6 for the whole PDG range of m a . Ac- 
cording to the scaling of Table [Til this supports the view 
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FIG. 11: (Color online) The meson constants f„ as functions 
of the o mass. 

that a is a qq meson, while /o(980) is a glueball. With 
this assignment the results of Fig. \TU\ emerge naturally. 
Certainly, the conclusion relies on the assumption of hav- 
ing just two states dominating the sum rules. Yet, it is 
appealing in its simplicity. In the following section we will 
see that if infinitely many states are considered this con- 
clusion would not necessarily follow when a further sum 
rule (the vanishing of a dimension-2 condensate) is im- 
posed. Note that in the two-state model the dimension-2 
object, C2, is necessarily positive, see Eq. 

Figure QT] shows the dependence of the constants f a 
and ft on m a . For comparison, we also plot the pion de- 
cay constant / w = 93 MeV. We note that f a is of the or- 
der /„• near m a ~ 500 MeV. As discussed in Appendix iBl 
/tt and f a need not be equal even if we postulate the 
scalar to be the chiral partner of the pion. 

Finally, we examine the last sum rule (1921) , concerning 
the energy density of the vacuum due to the gluon con- 
densate. According to Ref. [74|], we have (for three active 
flavors) 

* = -|t(-G 2 > = -(224i^ MeV) 4 . (94) 

We note that the two-state model favors the lower values 
of m a , up to 550 MeV, where the sum rule (p?2"j) is satisfied 
within the values (|94]l . 



VI. DIMENSION TWO CONDENSATES IN THE 
SCALAR SECTOR 

In the previous section we have explored phenomeno- 
logical consequences of saturating the sum rules in the 
large- A^ c limit with a finite number of states (two states) . 
Such an approach is justified for observables where the 
coupling constants to higher-mass states are sharply sup- 
pressed, and we are essentially left with the dominance 
of the low-mass states. This is, however, not always the 



14 



0.6 rr 

0.5 1 

> 0.4 1 
o 

o 

" 0.3 ; 

o.i 1 



QCD sum rules 



0.0 



0.4 



0.6 0.8 

liter [GeV] 



1.0 



1.2 



FIG. 12: (Color online) The fourth root of the minus energy 
density of the vacuum, (— e„) 1//4 in the two-state model as a 
function of m CT . The horizontal line and the band indicate the 
estimate with uncertainties of Eq. (|94[) f74j]. 



case. For instance, as discussed in Sect. IIIIB1 proper 
matching of the two-point functions to QCD necessar- 
ily requires infinitely many intermediate states present 
in the correlator. Only that way the leading logarithmic 
behavior ([2"2")) can be reproduced. Similar studies have 
been carried out in the past for other correlators in the 
Regge approach |75l485| . Here we focus on the Ilee cor- 
relator defined in Eq. (TTB| . 

In QCD, keeping the leading terms in a, this object 
has the explicit twist expansion of the form |5CH54| | 



n e e 



a 2 (N 2 



16tt 



32tt 
— G 2 ' 

7T 



C 2 Q 2 



(95) 



(note a factor of 16 difference between our definition (jTHJ) 
and Ref. [54], carrying over to the definition of C2). 

The dimension-2 gluon condensate, originally pro- 
posed by Celenza and Shakin [861 ] . appears as an elu- 
sive gauge-invariant non-perturbative and non-local op- 
erator which generates the lowest l/Q 2 power correc- 
tions in the twist expansion. Its dynamical origin re- 
mains unclear [§3, IH| , desp ite certain evidence provided 
by the instanton model [80], phenomenological QCD sum 
rules reanalyses [§0[ , phenomenological studies of the r 
decay data [9l|, further theoretical considerations (92| — 
|95| . quark- model calculations [H, [53], or lattice simu- 
lations [H, [9t| and their relevan ce fo r the confmement- 
deconfmement phase transition jl00| ] . Phenomenologi- 
cally, upper bounds for the dimension- 2 operato rs coming 
from the e + e~ data were analyzed in [I0ll4l03j . yielding 
16C* 2 < 0.1 GeV 2 - (0.3 GeV) 2 . The possible appear- 
ance of the dimension- 2 con densates in Regge-like models 
can be traced in Ref. [1041 ]; it was also briefly discussed 
in Ref. [8lj | . A more quantitative analysis was carried out 
by us in Refs. (8~ 



Very recently, the presence of the dimension-2 objects 
has be en l inked to the truncation of the full perturbativc 
series [l05j |. indicating duality between the resummed 
perturbation theory and the presence of the lowest power 
corrections. This is supported by the correlation between 
perturbative and non-perturbative contributions found in 
Ref. [l06| . 

As we have already discussed in Sect. IIII B[ in the 
6-6 correlator the dimension-2 object, Ci = ^2, n f 2 , 
appears naturally. Clearly, if there is a finite num- 
ber of states, the sum is necessarily positive and hence 
non-vanishing. This is in an open contradiction to 
the Operator-Product-Expansion philosophy, where only 
dimension-4 and higher local gauge-invariant operators 
are allowed. Therefore, the only scenario where the 
dimension-2 object can possibly vanish is by assuming an 
infinite number of resonances, such that the sum diverges 
requiring regularization. In this case the asymptotic con- 
dition, Eq. (|!Z3"]) . holds. As usual in the Regge models, we 
assume f n to be independent on n for all values of n, not 
only asymptotically, as requested by Eq. ©. 

The basic object of our study is the Regge sum (in this 
section M denotes the lowest glueball mass) 



n 



ee 



E 



4 i2 



Q f, 



14 fl 



J n 



ml + ' 



^ M 2 

n=0 



Q 2 ' 



(96) 



which can be evaluated using the formula (see Ap- 
pendix [AJ 



00 r 
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n=0 
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M 2 + an 

M 2 ~ 



Q 
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M 2 + Q 2 



(97) 



where 



T'(z)/T(z) 
log z 



1 



1 



2z I2z 2 I20z 4 



+ 0(z" 6 ),(98) 



is the polygamma function. 

We note that the infinite Regge sum generates properly 
the Q A log Q 2 term. We assume that we have one family 
of states contributing to the Regge sum at the leading- N c 
order. The matching with the leading term of the OPE 
expansion gives immediately 



Co 



11 

a 



a 2 (N 2 



32tt 4 



(99) 



We read off that / ~ N C) in accordance with the fact 
that the glueballs provide the saturation at the leading- 
N c level. Explicitly, 



/ = 



liy^v 2 

12\/27r 2 : 



(100) 



which for a = 0.2, a = 1.31GeV , and N c = 3 gives 
/ ~ 135 MeV, a reasonable value. Note that taking the 
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full slope a = 1.31 GeV is equivalent to taking every sec- 
ond state in the scalar-isoscalar trajectory, hence we are 
implicitly assuming that every second state on the tra- 
jectory is a glueball, alternating with qq scalar-isoscalar 
mesons. This alternation is suggested by the two parallel 
trajectories, Eq. (0. 

For the case of the dimension-2 condensate the £- 
function regularization [84[ (see Appendix [AJ applied on 



the spectrum M 2 



M 2 yields 



C 2 



(101) 



The condition for a vanishing dimension-2 condensate is 
thus equivalent to requesting that the mass of lowest state 
be given by the equality 



9 1 

- 2°» 



(102) 



which for a = 1.31(12) GeV from Eq. © means M = 
810(40) MeV. The value is not far from the mass of 
/o(980), it is also in the PDG range for m a (cf. Fig.©. 6 
Thus, on the basis of a simplest Regge model considered 
here it is not possible to sort out which state is a glue- 
ball, and which one is a qq meson, if the dimension-2 
condensate were indeed to vanish. 

The condition for the vanishing C2 is identical to 
that for the vector or axial flavor channels with Ciy = 
T, n fn,v and C 2,A = E n /n,A) respectively (Hf. Note 
that in that case the first Weinberg sum rule requires 
a non- vanishing dimension-2 object / 2 = C 2 ,a ~ C 2 y ■ 
This shows that both dimension-2 operators cannot van- 
ish simultaneously, despite the fact that they cannot be 
represented by local and gauge invariant operators. 

One may extend the Regge model by modifying the 
parameters of the lowest states on the trajectory, for in- 
stance the / constants. It is known from similar studies 
that such modifications cause strong effects [HJ. Also, 
the 1/N C effects may play a role. Nevertheless, we see 
that the fine-tuning needed to cancel the dimension-2 
object leads to a quite natural condition for the lowest 
mass of the Regge family. 

For the dimension-4 condensate the matching of 
Eq. to Eq. (|9"5j) leads to the large- N c identity (see 
Appendix |A"| 



l21N 2 a a 2 



2 



M a , r - 

+ - - M' 

a 



(103) 



The condition for the gluon condensate to be positive is 
therefore a > (3+\/3)M 2 or a < (3- V3)M 2 . Combining 



As described in 13211 , the shift from the pole resonance to the 
Breit-Wigner resonance is 0(Nc 2 ) suppressed, suggesting that 
the leading N c values lies in between. 
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FIG. 13: (Color online) The gluon condensate, (^G 2 ), in the 
one-family Regge model plotted as a function of the mass 
of the lowest state in the glueball trajectory, M. We use 
a = 1.31(12) GeV and a = 0.2(1). The band indicates 
the uncertainty in these quantities. The horizontal dashed 
line and the band indicate the estimate of Ref. [74|], equal to 
0.009(7) GeV 4 . 



Eq. fl03"l) and (fT00|) yields 



a 2 _ uNl (a 2 - 6aM 2 + 6M 4 ) 



(104) 



-G 2 ) = 0.022(3) GeV [10 



-G 2 ) = -0.005(1) GeV [10 



This quantity is plotted against the mass M in Fig. [13] 
We note that we have two windows for the values of M, 
above 1.3 GeV and below 0.45 GeV, which match the 
predictions of the considered Regge model to the QCD 
estimate of Ref. [z3 giving (fG 2 ) = 0.009(7) GeV 4 . Al- 
ternative determinations yield about twice larger values, 

or even negative values, 
. In our model the gluon 
condensate is related to the lighest glueball mass M ac- 
corging to Eq. (|103p . Modifications of the Regge model 
along the lines of Ref. [821 l84j may modify or lift these 
constraints, making it possible for the approach to be 
fully consistent with the requirements of the Operator 
Product Expansion. 



VII. CONCLUSIONS 

We have presented a wide-scope analysis of the cr-state 
in the large- N c Regge framework. Our main points are 
the following: 

1. The a state fits very well to the scalar-isoscalar 
Regge trajectory including all the known states 
from the PDG review. From the viewpoint of the 
spectrum, there seems to be no particular distinc- 
tion among the scalar-isoscalar states. The glueball 
or meson nature has to be sought based on the N c 
scaling of the couplings of the ++ states to the 
vacuum or to hadrons. 
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2. The mass of the lowest ++ state (the a) can also 
be extracted from the analysis of gravitational form 
factors of the pion and nucleon, available from the 
lattice calculations. Extracted masses are com- 
patible with the benchmark determinations within 
the expected accuracy of the large- N c framework, 
showing the consistency of the approach and the 
proximity of the large- N c limit to the real world in 
the scalar-isoscalar sector. 

3. The constraints of chiral symmetry and the high- 
energy behavior of the energy-momentum correla- 
tor imposes a set of sum rules involving various cou- 
pling constants. Saturation of these sum rules with 
just two states suggest that /o(600) is a meson and 
/o(980) is a glueball. Such an assignment, specific 
to the simple two-state model, results in natural 
iV c -scaling of the ratios of the coupling constants 
of these states. 

4. Correlation functions of the energy-momentum 
tensor are saturated by (infinitely many) scalar- 
isoscalar states. The analysis of the Operator 
Product Expansion shows that there appears the 
dimension-2 condensate, whose vanishing requires 
a fine tuning of parameters. Then, in the simplest 
Regge model the resulting mass of the lowest glue- 
ball state is ~ 800 MeV. 
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Junta de Andaluci'a (grant FQM225-05), the EU Inte- 
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Appendix A: The ^-regularization 



Appendix B: The linear Lagrangean for one scalar 



For only one scalar-isoscalar state, Lagrangean ([35 
with condition (|37l) reads 



A 
4/2 



(Bl) 



Introducing the Cartesian decomposition 



u 

fa 



<$>u 



(B2) 



one gets an exten sion of the linear sigma model with the 
so-called A-term [l09j |. allowed in the linear realization 
of the chiral symmetry: 



V(a 2 



where in the present case 



(B3) 
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Ja_ _ x 

J 7T 



(B4) 



Note that a and tt are chiral partners. 

Imposing the condition f a = / w removes the A term. 
Note, however, that the identity f n = f a is not a conse- 
quence of any symmetry or physical requirement. With 
this condition the Lagrangean reads 



In the large N c limit only tree-level diagrams survive, 
however, with infinitely many resonances. Thus, the infi- 
nite sums may need to be regularized. In the dimensional 
regularization the coupling of the resonance to the cur- 
rent acquires an additional dimension f n — > f n /J> £ with 
e = d — 4. By choosing as a natural scale u — m n , one 
gets f 2 — > f^m 2 ^, which means 



fn m 



2 k 
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2 m 2 J. 



(Al) 
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J 7T 



(B5) 



Note the resemblance to the bos onic part of the dilated 
chiral quark model of Ref. [llfj]. The main difference 
between potential (|B5|) and the standard Mexican hat 
potential, 



This is completely equivalent to formally expanding the 
two-point correlator at large Q 2 and reinterpreting the 
coefficients within the zeta-function regularization (84* 
In our case the mass formula is = an 
one has to consider the sums of the form 



l 0i 



(n- 



\2k _ 



lim > (n 



= C("2« 



(A2) 



which is the generalized Riemann zeta function. Then, 
we arrive at the explicit formulas 



V(a,if) 



^ 2 f 2 \ 2 
fir) i 



(B6) 



is the scale behaviour of the latter, a different vacuum 
energy density (factor of two difference), and modified 
couplings containing more than one scalar. The ctitit ver- 
tex is the same on shell, thus yields the same decay width 
of the a into two pions. 
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